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Abstract 

We examine W pair production in the Noncommutative Standard Model constructed with the 
Seiberg-Witten map. Consideration of partial wave unitarity in the reactions WW — > WW and 
e"*"e~ WW shows that the latter process is more sensitive and that tree-level unitarity is 
violated when scattering energies are of order a TeV and the noncommutative scale is below about 
a TeV. We find that WW production at the LHC is not sensitive to scales above the unitarity 
bounds. WW production in e^e~ annihilation, however, provides a good probe of such effects with 
noncommutative scales below 300-400 GeV being excluded at LEP-II, and the ILC being sensitive 
to scales up to 10-20 TeV. In addition, we find that the ability to measure the helicity states of 
the final state W bosons at the ILC provides a diagnostic tool to determine and disentangle the 
different possible noncommutative contributions. 
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I. INTRODUCTION 



With motivation coming from connections to string theory and possible relevance to 
quantum gravity, the possibility of noncommutative spacetime has received a great deal of 
attention in many different contexts in the last decade. The idea that spacetime coordinates 
might not commute, however, dates to Heisenberg [1], who originally hoped that quantum 
field theory on a noncommutative spacetime (NCQFT) could be free of the ultraviolet di- 
vergences that typically plague QFT. This motivated Snyder in 1947 to furnish the first 
explicit construction of a noncommutative spacetime [3]- The idea that NCQFT could cut 
off divergences dwindled since the success of renormalization, though it has received some 
attention in recent years 3[[. The theory and phenomenology of NC quantum mechanics has 
been well studied Q, Is], [g, [t]. Connes and his collaborators have sought to derive the Stan- 
dard Model gauge groups and representations as a unique consequence of NC geometry 



There is a strong re^ 
by many authors [ol. 



ationship between NCQFT and string theory which has been studied 
1^. 



In particular, interest in noncommutative extensions to the Standard Model increased 
when Seiberg and Witten [ll] described how NC gauge theory could arise as the low-energy 
limit of string theory with a background field. This leads to a commutation relation between 
spacetime operators, 

[x^'^x•']=^e^"', (I.l) 

where here, and throughout the rest of the paper, hatted quantities are in the algebra of 
operators generated by the non-commuting coordinates, and 6'^'^ is a constant, antisymmetric 
matrix. Only "space-space" noncommutativity {6oi = 0) can arise as a consistent limit of 
string theory in Seiberg and Witten's construction. One can, however, try to formulate a 
"space-time" (6'ot 7^ 0) NCQFT nonetheless; though this leads to potential problems as we 
will discuss below. 

There are several approaches to constructing a NCQFT. Starting with the canonical 
commutation relation, one can use the analogue of Weyl quantization to associate a function 
/(x) to an operator / by 



where / is the Fourier transform of /. One can then show [10| that the product of two 



/ = / 7i^/(^)e^'''''' (1-2) 
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operators fg = f-*^g, with 



fix) * g{x) ^ fix) expi-e^'' d , d M^) 



(1.3) 



being known as the Moyal star product. The commutator hkewise has a noncommutative 
version 

[f{x) t gix)] = fix) -kgix) - gix) -k fix) , (1.4) 

known as the Moyal bracket. 

Using this, one can estabhsh a map between any NC Lagrangian buih from operators 
(hatted fields) and an ordinary Lagrangian built from conventional fields. The NC La- 
grangian is formulated with star products replacing ordinary products and Moyal brackets 
substituting for commutators. It is easy to show that J d^xf -k g = j d^xfg, so it is clear 
that this procedure leaves kinetic and mass terms unchanged from the commutative case, 
but introduces momentum- dependent phase factors in the interaction terms. In a gauge the- 
ory, the Moyal bracket also leads to new interaction terms, such as three- and four-photon 
vertices in NC quantum electrodynamics. The primary impediment to formulating a non- 
commutative Standard Model (NCSM) in this way is that SUiN) groups do not close under 
star multiplication. UiN) groups do, however, and an NCSM has been constructed using 
a f/(3) X f/(2) X U(l) gauge group with some new fields added to implement symmetry 
breaking effectively 12(]. This NCSM, however, was shown to violate unitarity in gauge 



boson scattering at high energies 



13| 



In a different approach, Seiberg and Witten ll| showed how to solve NC and normal 
gauge equivalence conditions and derived expressions for noncommuting gauge fields in terms 
of normal gauge fields, order by order in 6^'^ . Called the Seiberg- Witten map (SWM), this 
construction is viable for SUiN) gauge groups. The SWM formulation of the NCSM thus 
uses the SM gauge group and particle content, but adds new interactions Q,Q,[l6l. This 
is the model we will study in this paper, and we will simply refer to it as the NCSM from 
now on. 

Because of the nonlocal nature of Eq. 11.31 the quantization of NC field theory is fraught 
with subtlety. Many authors have pointed out potential issues with unitarity, especially 



in the space-time NC case \Vn, ll 
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23| . On the other hand, with suitable 



modifications to the formulation of perturbation theory, others have claimed that unitary 



space-time NCQFTs can be constructed 
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271, 
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30|. We will not review 



;his discussion in any detail, but instead will take the point of view espoused by Calmet 



3l| , that the NCSM formulated using the SWM is a valid effective theory with a Hermitian 



Hamiltonian which does not intrinsically violate unitarity, but may violate partial-wave 
unitarity at some energy scale. ^ Inspired by the above study of unitarity in the U{3) x f/(2) x 
f/(l) model |13|, here, we will determine that scale by studying partial-wave unitarity in 
e~^e'~^W~^W~m the NCSM. We will find that the terms in this amplitude that grow with the 
center-of-mass energy do not completely cancel as in the SM. At a certain energy, therefore, 
partial-wave unitarity is indeed violated. We will then discuss the ensuing implications for 
collider phenomenology of this model. 

There has been a great deal of work on the phenomenology of the NCSM. Some very 
strong bounds have been placed on NC physics in general, but whether they apply to the 
NCSM formulated using the SWM has not, in every case, been established 



3j] . Defining 
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IIP 



/A^, where c^y is an antisymmetric matrix of 0(1) numbers and A is the scale 



where NC effects become relevant, the strongest bound on the NCSM comes from Lorentz 



violation 



36 



37 



38 



39 



Many possible 



40 



41 



42 



43 



phenomenological signatures of this model have been examined 



by Alboteanu, Ohl, and Riickl 
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5l[ |. Here, in a similar vein to work 



531], we study the phenomenology of the NCSM at LEP, 
the LHC, and the ILC. We concentrate on W~^W~ production, expecting that the terms in 
the amplitude that grow large at high energy and threaten unitarity will also lead to strong 
experimental signatures of the NCSM. We find that this is the case, and that by studying 
this process one can discover the NCSM and measure its parameters at future colliders. 



II. THE NONCOMMUTATIVE STANDARD MODEL 

In the Noncommutative Standard Model, the action is expanded in powers of the anti- 
symmetric noncommutativity tensor 6'^'^ using the Seiberg-Witten map. No new particles 
are present in the model, but there are new interactions and modifications to existing inter- 
actions. 

The Seiberg-Witten map relates each NC field to a nonlinear function of ordinary fields 

We note that the this model is in fact renormahzeable with a shghtly modified Lagrangian. 32] It would 
be interesting to see how that modification affects the phenomenology. See [s^, p. 42 for a more detailed 
discussion of this issue. 
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and 9^^ . The form of this function, for each field, is obtained by solving the so-called gauge 
equivalence condition, the requirement that an NC gauge transformation of NC fields induces 
an ordinary gauge transformation of the associated ordinary fields. For example, for an NC 
gauge theory with gauge field A, matter field "ij), and gauge parameter A, the SWM maps these 
NC objects to the corresponding functions A{A,6), 'ip^ip^A^O), and X{\,A,6), respectively. 
The NC gauge transformations have the same form as ordinary gauge transformations, with 
products replaced by star products as expected, 

^ = exp{ig\*)ijj = 'il) + igX-kil) + -^^A -k\i<il) + &{\^) , 
Af, A\^ = exp{igX'k)A^exp{-igX'k) + -exp{igX'k){df,exp{-igX'k)) (II. 1) 

= A^ + ig[X ^ A^] + i|^[A t [A t A,]] + d,X + ig[X t d,X] + @{X') . 

The gauge equivalence conditions then state that 

A{A,e)=A{A,9) , 

^'(^,A^)=V^(V^',^',^) (11.2) 
A'(A,yl,^) = A(A',A',^) . 

Solving these conditions to first order in 6^'^ yields 

A^{x) = A^{x) + ^eP^{A„{x), dpA^{x) + F,^{x)] + 0{e^) , 
^(x) = ^{x) + ie^M.(a;)9p^(x) + ^^'^[^.(x), A^(x)]^(x) + 0{d^) , (II.3) 

Z o 

A(x) = A(x) + ^^^'^{A.(x), d,X{x)} + &{e^) . 

The NCSM Lagrangian can then be formed, essentially, by starting with the SM Lagrangian 
and replacing all fields with their NC counterparts. Then, by using the SWM map given in 
the above expressions and Eq. II. 3t the Lagrangian can be expanded order- by-order in 9^^ . 

There are two caveats to this simple procedure. The first is that in order to solve the 
charge quantization problem that plagues NC gauge theories js^, one must introduce a 
separate NC gauge field for each U (1) eigenvalue. Since each of these is mapped to the same 
SM gauge field under the SWM, however, there are no new fields in the effective theory jl^ . 

The second is that, as can be seen by the anticommutator on the right-hand side of the 
first line of Eq. III.3t the NC gauge field A takes values in the enveloping algebra of the Lie 



algebra of the SM gauge group. This means that the SWM does not fully determine the 



gauge kinetic terms. The resulting freedom in the gauge sector is parametrized [4^ with the 
constants ki, K2, and K3. Since only K2 contributes to the interactions of the W boson, it is 
the only parameter relevant here. The choice = ^2 = = is called the minimal NCSM 
(mNCSM). Throughout this paper we will generally work with the non-minimal NCSM, 
considering various values of K2 and determining bounds on it and A from collider data. 
There are also bounds on the k parameters from consistency conditions on the theory 42 1. 

As mentioned above, the scale of noncommutativity A appears as 6'^,^ = c^,^/ h?^ where in 
this paper we will take Cp^p to be an antisymmetric matrix with each element being either 0, 
1, or -1. The space-time case Cqj 7^ 0, as stated earlier, has been shown to be problematic for 
unitarity, and we will see below that it is only this case that leads to dangerous unitarity- 
violating terms in WW scattering. We will also find that the space-space elements of 9^" 
give significant contributions to specific W polarizations in WW production at colliders. 

Here we present the Feynman rules relevant to our discussion as derived in the NCSM to 
&{9). In the diagrams displayed below, a vertex with a square refers to the @{9) contribution 
;o that Feynman rule. Most of the rules relevant to this analysis are derived in Melic et al. 
isl All gauge boson momenta are incoming. 

First, there are modifications to the SM fermion-gauge boson 3-point couplings. 




4^2 sin 9w 




^ Vp^V , (11.4) 
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where 6'/,^^ = 6'^^7p + O^p'y^ + Op/^'yu ■ 

Then there are modifications to the SM three-gauge-boson vertices, 



2 f iJ.vpi.P^ 

+ 2esm26wKww^Oixiyp{p,q,r) , (II.5) 



+ 2esm29wKwwz0^^p{p,q,r) , 

where 

f^upiP) =0t.uPp + d,,pPi^ 

+ 9p.v{0 ■ p)p - gvp{0 ■ p)m + 9pn{d ■ P)v , 
fpupiP, = 0^u{p - q)p + 0^p{q - r)^ + Oppir - p)^ 

-2gp,{e ■ r)p - 2g,p{e ■ p)^ - 2gpp{9 ■ q), , (II.6) 
Oi,up{p, q,r)=e^,„{p-r qp-q-r Pp) + {9 ■ p)p (q -r g„p- qpT^) 

- {6 ■ P)v {q -r 9pp- qpVp,) - {9 ■ p)p {q -r g^,y- q^,r^) 

+ px q{rpg^p-r„gpfj) 

-\- (cychc permutations of {p, g, r} and {/x, p} simultaneously) , 

Kwwj — — 5'^/^2/2c^s^, Kwwz — 5'^'«2/2c^, and p x q = p'^9pi,q''. Here, K2 is the parameter 
capturing the freedom in the gauge sector discussed above. 
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There is also a 4-point fermion-gauge boson interaction, 




= -^{o,.p{p' + K)ii - 75) + 2V^) ■ 



(II.7) 



Finally there is a modification to the 4-W vertex, which was not previously presented in 
the literature. 



X 

Note that all of the 0{O) contributions to the Feynman rules are opposite in phase from 
the SM contributions (i.e. one is real, while the other is imaginary). This fact determines 
the lowest order in 9 at which the NCSM corrections to coUider observables occur. To see 

why this is, we can examine the phases of various contributions to a scattering amplitude. 
For example, consider W^W^ W^W^ scattering with s-channel photon exchange. Up 
to an overall phase, the amplitude can be written as 

A = eMVs'M + ^VZ'){^)iVs'Z' + ^VZ')e;e; , (119) 

where e is a polarization vector for one of the external W bosons, and V, which represents 

the 'yW~^W~ Feynman rule, can taken to be real. Vsm is the Standard Model Feynman rule, 
while Vnc is the NC correction at first order in 9. Truncating the amplitude at first order 
in 9 yields 

A - ^eMiVs^MVs'M' + KVs^mVZ' + VrM^V^'Mel . (11.10) 

Defining = ef.e^V^^ yields V^* = e*e*V^''^, since V is real. The amplitude can now be 
written as 

A^--{VsM- V*sM + i{VsM ■ VjVc + V*sM ■ Vivc)). (11.11) 
s 

If the vectors are real, then V is real, in which case the squared amplitude is 

1^1' = ^{(VsM ■ Vsm)^ + (Vsm ■ Vnc + Vsm ■ Vnc)^), (11.12) 

8 



which has no terms that are first order in ^. In fact, so long as each has only an overall 
phase, that is, each Lorentz component of a given polarization four-vector has the same 
phase, then the phase can be factored out into an overall phase multiplying the whole 
amplitude and the argument still goes through. To obtain a first-order contribution, we 
require that at least one of the polarization vectors has a relative phase difference among its 
components. For W~^W~ W~^W~, this means that at least one of the W bosons must 
be transversely polarized. In the sum over all polarizations, however, 0{9) contributions to 
the amplitude must vanish. 

With fermion external states, however, every helicity choice has a spinor with relative 
phase differences among its components. In this case there are nonzero 0{O) contributions 
to the squared amplitude for any helicity choice and for the sum over helicities. 

III. PARTIAL WAVE UNITARITY 

It is expected that any perturbative model of new physics satisfy tree-level unitarity. 
It is well-known that any amplitude that exhibits azimuthal rotational invariance can be 
decomposed into partial waves ai according to 

oo 

A= 167r^a,(2Z + l)P,(cos0), (III.l) 

1=0 

where 9 is the scattering angle in the center-of-mass frame. Pi is the Ith Legendre polyno- 
mial, and ai is called the Ith partial wave amplitude. The orthogonality of the Legendre 
polynomials then gives 

1 

ai^ / A{cos9)Pi{cos9)dcos0. (III.2) 

327r J_i 

It is a consequence of probability conservation that 

Re{ai) < 1/2 , < Im(a/) < 1 , and |a,|^ < 1 . (III.3) 

In the Standard Model without a Higgs boson, the W'^W'^W^W' amplitude contains 
terms that grow linearly with s. At high enough scattering energy, these terms grow large 
enough to violate the above partial-wave unitarity bounds. The inclusion of Higgs exchange 
diagrams cancels these growing terms, and (for large Higgs mass uih) the leading terms that 
remain are proportional to m^. Partial- wave unitarity then gives an upper bound to the 
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Higgs mass tuh ^ 800 GeV. One can also check unitarity in other processes in the SM, such 
as e+e^ W~^W~. In this case the terms in the amphtude that go as s once again cancel, 
though they do so in the absence of diagrams involving the Higgs, as they should, since the 
Higgs diagrams vanish in the limit of small electron mass. 

If, in a model of new physics, this cancellation of terms that grow with scattering energy 
does not occur, unitarity will be violated at some scale. Above this scale, the model is 
invalid or incomplete. As we will see, this is the case in the NCSM, though since a partial 
cancellation occurs for W~^W~^W^W~ , a stronger effect occurs in e~^e~^W~^W~ . 

The partial wave unitarity analysis in the NCSM is complicated by the fact that the 6^1, 
tensor gives preferred directions in space. This breaks the azimuthal rotational invariance. 
The simple partial-wave analysis discussed above is thus not necessarily valid. The appropri- 
ate form for the partial-wave unitarity bounds for amplitudes without azimuthal symmetry 



was given by Chaichian, Montonen, and Tureanu [SOj. These bounds apply to a partial- wave 
expansion of the amplitude in all of the independent angular variables. Evaluating these 
bounds requires calculating the amplitudes with an arbitrary crossing angle, which is beyond 
the scope of this paper. 

To get an estimate of the unitarity constraints, we will instead remove the (p dependence 
from the amplitude by fixing at a value that maximizes the NC contribution and use the 
bounds for azimuthally symmetric systems given in Eq. IIII.3I While not strictly accurate, 
if we fix (f) so as to maximize the NC contribution, this procedure should provide a "worst- 
case scenario." If the unitarity bounds are weak in this case we can conclude that a more 
thorough analysis will not give significantly stronger limits. 

We will first examine W~^W~ ^W~^W~ , the process that is most dangerous in the SM. 
In the NCSM, there are twenty diagrams contributing to this process at tree level. There 
are six s-channel gauge boson exchange diagrams. 






(111.4) 
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six t-channel gauge boson exchange diagrams, 




(III.7) 



and, finally, two contact interaction diagrams. 





(III.8) 



Here again, a vertex with a square refers to the 0(6*) NCSM contribution. 

We will consider the case in which all four W bosons are longitudinally polarized, as 
(like in the SM) this gives the worst high-energy behavior. We also assume a coordinate 
system where the beam line is aligned along the z- or 3-axis. We note that in fact, if the 
orientation of 6 were fixed with respect to, say, the cosmic microwave background, then the 
earth's motion would continually alter the orientation of the beamline with respect to 6. It 
is a straightforward though involved procedure to account for this (it was carried out for a 



study of noncommutative QED at LEP 



55|), and will be neglected here. 
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For the SM contribution, as mentioned above, the leading terms are independent of s, 
and in the hmit mn ^ simphfy to 

2 2 

Msm{W,^W,- iyo'W^o")leading terms in s = , (111.9) 

Where the subscript "0" indicates longitudinal polarization. For the NCSM contribution, the 
leading terms for each diagram go as s^, but these cancel in the sum over diagrams. There 
is an incomplete cancellation among the terms that go as s, however, and what remains is 



(cos^ + l)(cos^-3) 



sin^^(cos0coi + sin0co2) + (cos^^ — l)co3)s . (111.10) 



164A2 cos^-l 
Setting (p = tt/2 and using Eq. 1111.21 we get 

«o(»VWo- - W-W,-) = + .55^(32c,3 - 21.e„). . (III.U) 

to which we can now apply the bounds of Eq. 1111. 3[ Because Mnc is pure imaginary, the 
bounds on the magnitude and on the imaginary part of ao are the relevant ones, and they 
give numerically equivalent results. We choose Co2 = 1 and c^i, = otherwise to get the 
strongest bound on A, shown in Fig. [H Only very small values (compared to ^/s) of A are 
excluded. 

Since partial-wave unitarity in W^W~ ^W^W~ does not give a meaningful bound on 
A, we now turn to the process e~^e~^W~^W~ . There are ten diagrams that contribute to 
this process at tree level in the NCSM. First, there are s-channel gauge boson exchange 
diagrams. 






(111.12) 

Second, there are neutrino exchange diagrams. Note that for up-type fermions these will 
be ti-channel diagrams, while we show the t-channel diagrams appropriate for down-type 
fermions. 






(111.13) 
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EcM (GeV) 



FIG. 1: The bound on the noncommutativity scale A as a function of i/i from partial-wave 
unitarity violation in WW scattering. In the shaded area, unitarity is violated. 

Finally, there is the contact term. 




(III. 14) 

We will examine the helicity choice e^e^ WqW^^', which has the worst high-energy 
behavior. Adding the diagrams, we find, as mentioned above, that the terms which grow 
with s cancel in the SM, and the remaining constant term is numerically negligible. Un- 
like for W^W~^W~^W~ , however, in this process the NCSM terms that go as do not 
completely cancel, and we are left with 

■^(e^efl Wq iyo^)leadmg terms in s = o 2 A 2 ^^^03 ^ ^ ' (III.15) 

We once again employ our procedure regarding 0, in this case we set = and use 
Im(ao) < 1/2. Because of the dependence, we find strong bounds on A as s increases, 
as shown in Fig. O This unitarity limit should be compared to the collider search reach in 
NCSM parameter space, which we discuss below. 
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Unitarity bound for co3=l from e^eJ^V^o Wq amplitude 




EcM (TeV) 

FIG. 2: The bound on the noncommutativity scale A as a function of y/s from partial-wave 
unitarity violation in e^e~—>-W~^W~. In the shaded area, unitarity is violated. 

IV. COLLIDER PHENOMENOLOGY 

Since unitarity does not exclude the NCSM as a valid or complete theory at LHC and ILC 
energies, it is natural to investigate its phenomenology at these colliders. Here we examine 
its effects in WW production. We will show that this process gives unique signatures that 
are easily discernible for large regions of parameter space, and that by measuring the W 
polarization one can disentangle and measure the different parameters of the model. 

As mentioned above, our calculations of collider observables assume that the 3 direction 
[i.e. the direction corresponding to = 3 or z/ = 3 in 6^^,^) is aligned with the collider's 
beam axis, and we do not undertake the procedure fssi of taking into account the earth's 
motion with respect to the frame in which 6'^,^ is fixed. 

A. LEP phenomenology: e+e~ — > W'^W~ 

The process e~^e~ W~^W~, besides giving unitarity constraints on the NCSM as de- 
tailed above, also provides observable signatures at high-energy e+e~ colliders. The NCSM 
predicts deviations in the e~^e~ — >■ W~^W~ differential cross section that, for certain regions 
of parameter space, would have already been observable at the high-energy running of the 
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Large Electron-Positron Collider. 

Each W boson can decay either hadronically or leptonically. There are then three distinct 
final states for e^e^ W~^W~, hadronic semileptonic, and leptonic, which occur 45%, 45%, 
and 10% of the time, respectively. At lepton colliders, including LEP and ILC, each of these 
final states is easily observable with high efficiency, and in the hadronic and semileptonic 
cases the full event can be reconstructed. For the LEP analysis, we can therefore simply 
calculate the e~^e^ —>■ W^W~ differential cross section and multiply by 90% to exclude 
the fully leptonic events where reconstruction is impossible and therefore the W~ direction 
cannot be determined. 

We calculate the amplitude for this process using the diagrams displayed above. After 
squaring the amplitude and summing over helicities, we find that for the 0(6') correction to 
the differential cross section, the dependence on the noncommutativity tensor 9^^ and on 
the azimuthal angle takes the simple form 

6da /d cos 9 d(f) (X —(cqi sine/) — co2 cos 0) . (IV. 1) 

A^ 

Thus LEP-II is only sensitive to space-time noncommutativity. Also, any observables that 
are integrated over the full range in (p will not be sensitive to the NC contributions, since 
COS0 and sin0 both integrate to zero. Also note that the leading term in the amplitude 
given in Eq. IIII.15I does not contribute to the squared amplitude at leading order in 9. 

To determine the region of parameter space excluded by LEP-II, we take the double 
differential cross section 6da/dcos9 d(f) binned in 20 x 20 equally sized bins. We calculate 
the fo^^ this distribution with respect to the SM, assuming a total integrated luminosity 
of 700 pb^^ and taking into account statistical error plus a 0.1% luminosity uncertainty. 
Choosing cqi = co2 = 1, we vary A and K2 to find the exclusion contours. 

LEP-II operated at a number of center-of-mass energies ranging from 130 GeV to 
209 GeV. Running at multiple energies is essential to rule out acceptance effects as an ex- 
planation for any putative azimuthal dependence. Azimuthal dependence from new physics 
should change (typically, increase) with increasing energy, while azimuthal dependence from 
acceptance effects should not change with energy. To simplify our analysis, however, we just 
take all of the integrated luminosity at Ecm = 200 GeV, which is close to the luminosity- 
weighted average scattering energy. 

Fig. [3] shows the resulting exclusion region from W pair production at LEP-II. The reach 
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in A is above the center-of-mass energy for all values of k,2- 

LEP exclusion region for coi=co2=l 

4 
2 

^ 
-2 
-4 

100 200 300 400 500 600 

A (GeV) 

FIG. 3: From left to right, The five-cr, 95%, and 68% confidence level exclusion contours 
from the process e'^e~ —>■ W^W~ for the NCSM at LEP-II. We also display the bound 
obtained previously from partial-wave unitarity. 

A typical example of the signal that would be observed in the excluded region is given 
in Fig. m Here, as throughout the paper, we plot the differential event rate dN/ d cos 9 d(j), 
normalized such that the value in each bin is the actual number of events that would fall in 
that bin. Shown is the distribution for the 0.8 < cos^ < 0.9 bin. The error bars depict 
the same errors used in the analysis-statistical error plus a 0.1% luminosity uncertainty. 
For A = 150 GeV and ^2 = 1, the distribution deviates from the SM well beyond the 
errors, and the characteristic sinusoidal shape is visible. 

B. LHC phenomenology: pp^W~^W~ 

At the LHC, W pairs are copiously produced (in both the SM and NCSM) via the process 
qq W~^W~ , with q = u,d giving the dominant contribution. The NCSM diagrams for 
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FIG. 4: A comparison of LEP-II data for the SM (dotted) vs. the NCSM (sohd). Shown is 
the the 0.8 < cos^ < 0.9 bin of the double differential cross section dN /d cos 6 d(j) for 
e~^e~^W~^W~ . The error bars are statistical plus certain systematics (see text). 

these processes are the same as were given above for e+e^ — » W'^W' , except that for uu 
initial states, the t-channel diagram is instead a w-channel diagram. The squared amplitude, 
once summed over helicities, exhibits the same 6^^, and dependence as in e^e~ W~^W~, 
given above. 

The large QCD backgrounds at the LHC swamp the hadronic final state, and, as at 
LEP, for the leptonic final state the event cannot be reconstructed. Because the relative 
momentum of the center-of-mass frame is unknown, the full event reconstruction in the 
semileptonic case is also not possible at the LHC (unlike at an e^e~collider). If one makes 
the approximation that the leptonically-decaying W is always on shell, however, then there 
is enough information to reconstruct the event. Quantifying the validity of this approxi- 
mation would require calculating the full 4-fermion production amplitudes and using these 
to generate events. This is beyond the scope of this paper. Here, in order to determine 
the LHC sensitivity to NCSM parameter space, we calculate the production rate of stable, 
on-shell W bosons and assume that reconstruction is possible in the semileptonic 45% of 
events. Because the W width is relatively small, this should give a fairly good approximation 
to the true search reach. 



To obtain the differential cross section, we use the Mathematica implementation 5^ of the 
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CTEQ5 parton distribution functions [57| • Because, for a hadron collider, the longitudinal- 
boost-invariant WW invariant mass mww is a more convenient variable, we will consider 
the WW invariant mass distribution dN/ dmwy/dip. This also takes advantage of the spread 
of scattering energies available at the LHC. The NCSM effects become easier to observe at 
higher invariant masses where the the SM event rate falls off. We choose 20 equal-sized bins 
spanning from threshold, my/w = 2mvi/, up to mww = 5.5 TeV. 

Once again we calculate the for the NCSM distribution relative to its SM counterpart. 
As we did for LEP, we include statistical error and a luminosity uncertainty, which we take 
to be 5% for the LHC. Taking coi = co2 = 1, we use the to determine the LHC discovery 
reach in the A.-K2 plane. We assume an integrated luminosity of 100 fb^^. The result is 
shown in Fig. O 




FIG. 5: The search reach for the NCSM from pp W+W' at the LHC for 100 fb"^ 
integrated luminosity. Superimposed are the perturbative unitarity bounds for specific 

parton scattering energies. 

Superimposed on these exclusion contours are the unitarity constraints on A for various 
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values of ^/s, obtained from Fig. [2l For example, in order for the model to satisfy tree-level 
unitarity at scattering energies of Ecu = 1 TeV, A must be greater than about 200 GeV. 
The LHC search reach does not extend very far into the region of parameter space where 
the NCSM is unitary at typical LHC partonic scattering energies, so we conclude that the 
process pp — > WW is not an optimal tool to probe the NCSM. 

A typical example of the signal that would be observed in the discovery region is dis- 
played in Fig. [61 Shown is the distribution in the 1.4 < mww ^ 1-6 TeV bin of the 
double differential event rate dN/ dmy/wdip. The error bars depict the statistical errors plus 
luminosity uncertainty used in the x^- For a low enough value of A, the sinusoidal deviation 
from the fiat SM prediction is observable. 

SM vs. NCSM event rate for LHC, 100 fb"' 
A=600 GeV, coi=co2=l, K2=3, mwwG[l-4,1.6] TeV 



100 
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1 2 3 4 5 6 

FIG. 6: A comparison of LHC expectations for the SM (dotted) and the NCSM (solid). 
Shown is the 1.4 < m^w ^1-6 TeV bin of the double differential event rate 
dN / dniy/wdct). The error bars are statistical plus certain systematics (see text). 

The signal should also be visible in the invariant mass distribution, so long as we choose 
a certain range in 0. From Fig. [6] we can see that integrating over the range from 1 ^ ^ 4 
should maximize the excess of the NCSM over the SM. Figure [7] shows the invariant mass 
distribution with the integral taken only over this range. A systematic excess of the NCSM 
over the SM is apparent at high values of the invariant mass. 



19 




FIG. 7: A comparison of LHC expectations for the SM (solid) and the NCSM (data 
points) for the differential event rate dN/dmy/w- 

C. ILC phenomenology: e'^e^ ^W^W~ 

Our analysis for the ILC is similar to that for LEP except that the center-of-mass energy 
is higher (500 GeV and 1 TeV), the luminosity we consider is greater (500 fb~^), and the 
luminosity uncertainty is smaller (10~^). Also, the ILC will have polarized beam capability, 
so in addition to the differential cross section da /d cos 6 d(f), we can consider the differential 
left-right asymmetry dAin/ d cos 6 d(f), where the left-right asymmetry Alh is defined as the 
asymmetry between cross sections with fully polarized electron and positron beams, that is 

^ _ cr(e^e+ ^ W^W~) - a(e^e+ ^ W^W~) 

Of course, at a real collider, the beam polarization is imperfect, but is straightforward to 
relate the asymmetry measured with partially polarized beams to the "theoretical" value of 
Alr defined above. The beam polarizations and the accuracy to which they are known of 
course affect the errors on the measurement of Alr and thus the bounds obtainable from 
it. We assume that the electron and positron beams have polarizations Pg- = 0.9 and 
Pe+ = 0.6, respectively, and that the polarization uncertainty is AP/P = 0.25%. It should 
be noted that Alr shows the same dependence on 9^i, and as the unpolarized differential 
cross section. As we will see in the next section, distinguishing contributions from the other 
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components of 6 requires measuring the W polarization. 

Using these two observables, we again perform a analysis and determine the exclusion 
contours for the ILC. Figure [H] shows the discovery reach for the ILC assuming 500 fb~^ of 
operation at ^/s = 500 GeV, taking cqi = Co2 = 1- Shown are contours both with and 
without the inclusion of dAm/ d cos 6 dcj), along with the unitarity bound for ^/s = 500 GeV. 
Unlike at the LHC, at the ILC the search reach extends well beyond the unitarity bound. 
Also evident is a dip in the no-ALpt exclusion contours around /€2 ~ 2 where the unpolarized 
observables are relatively insensitive to the NC effects. This is due to a partial cancellation 
in the amplitude between the minimal and non-minimal NCSM contributions. As can be 
seen from the other set of contours, the addition of A^ji as an observable untangles this 
cancellation and improves the search reach greatly. 

500 GeV ILC exclusion region 
coi=co2=L Xint=500 fb"^ 

4 
2 

5^ 
-2 
-4 

0123456789 10 

A (GeV) 

FIG. 8: From left to right, The five-cr, 95%, and 68% confidence level discovery reach for 
the NCSM from e+e-^W+W-at the 500 GeV ILC, with and without Alr. 

Figure [9] shows a typical example of the unpolarized signal in the discovery region. We 
display the -0.9 < cos^ < -0.8 bin of dN/dcos9d(t) for A = 1 TeV and = 1- The 
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sinusoidal deviation from the flat SM distribution is apparent. 




FIG. 9: A comparison of 500 GeV ILC expectations for the SM (dotted) and the NCSM 
(solid). Shown is the the —0.9 < cos^ < —0.8 bin of the double differential cross section 
dN /d cos 6 d(j) for the process e^e~ ^W^W~ . The error bars are statistical plus certain 

systematics (see text). 

In Fig. [TOl the A^^ signal is presented for a point in parameter space (A = 3 TeV, 
^2 = 2.4) where the unpolarized differential cross section has negligible deviation from the 
SM. The differential A^ji shows, however, a clear sinusoidal deviation characteristic of the 
NCSM. 

We also consider the process e^e~^W~^W~at a y/s = 1 TeV ILC with an integrated 
luminosity of 500 fb^^. As can be seen in Figs. [TT], [121 and [13], the phenomenology is 
qualitatively similar to the ^/s = 500 GeV case but the search reach is significantly higher. 
In both cases, the ILC, especially when utilizing polarized beams, is sensitive to values of A 
many times greater than ^/s. 



D. W polarization at the ILC 

All the observables discussed above have the same dependence on 9fj_^ and (j), that is, 
oc 6^01 sin — 6*02 cos 0. This is true of any observable that does not distinguish between 
final state polarizations. At the ILC, however, W polarization can be determined with high 
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FIG. 10: A comparison of 500 GeV ILC Alr data for the SM (dotted) and the NCSM 
(sohd) for e~^e~^W~^W~ . Shown is the the —0.2 < cos 6* < —0.1 bin of the double 
differential left-right asymmetry d A ^ ^ / d cos 6 d(j). The error bars are statistical plus certain 

systematics (see text). 



efficiency [58|. By comparing the production of different combinations of W polarizations, 
we can disentangle the contributions from different elements of Of^^, and from K2- In this way, 
we can measure these different parameters at the ILC. 

In this section, we will compare differential event rates for W^W~ production at a \/s = 
1 TeV ILC with a 500 fb~^ integrated luminosity for different and W~ helicities. To 
compare different NCSM contributions, each NC distribution we present corresponds to one 
component of c^^^ = 1 with all others are set to zero. The analytical expressions for the 
amplitude squared for the various helicity choices are given in the appendices. 



1. dN/d(j) 

First we examine the differential event rate dN/ d(f), integrating over the full range in cos 6. 
Figure [TH compares W^W~ production at a 1 TeV ILC for the SM and the NCSM for both 
{R,0} and {0,L} helicity choices (for {l^"*",!^"}), and for both ^2 = and ^2 = 1. Here, 
refers to longitudinal polarization, and R and L to right and left circular or transverse 
polarizations, sometimes denoted by + and -. As is clear from the plot, with A = 3 TeV, 
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1 TeV ILC exclusion region 

C01=C02=1, -Cint=500 fb"^ 




5 10 15 20 

A (TeV) 



FIG. 11: Same as Fig. El except for = 1 TeV. 




FIG. 12: Same as Fig. [9l except for = 1 TeV, A = 2 TeV, and -0.9 < cos^ < -0.8. 



24 




FIG. 13: Same as Fig. ^ except for = 1 TeV, A = 5 TeV, and -0.7 <cos9 < -0.6. 

the Ci3 = 1 and C23 = 1 distributions are easily distinguishable from each other and from 
the SM. The K2 = 1 NCSM distributions are also clearly distinguishable from the K2 = 
distributions. 

As can also be seen, the NCSM contribution from C13 = 1 is opposite in phase for the 
{R,0} and {0,L} helicity choices. When summing over these two helicity choices, the C13 = 1 
contribution vanishes. The same is true for the C23 = 1 contribution. In order to determine 
the values of C13 and C23, therefore, it is essential to measure the W polarizations. 

Different information can be extracted if we look at the dN/ dcf) distribution in a particular 
cos 9 bin, instead of integrating over all values of cos 9 as we did in Fig. [T3] above. Figure [T5] 
shows the W^j|iyQ~ distribution for only the —1 < cos9 < —0.9 bin. Here, we see that the 
C12 contribution is distinguishable from the SM and other mNCSM (^2 = 0) contributions 
(which are negligibly different from the SM in this case). This distribution is also sensitive to 
the value of k,2] the K2 = 0, distribution is significantly different from the K2 = 1 distribution. 
Interestingly, this contribution is flat in rather than having the sinusoidal shape that has 
characterized the NCSM distributions we have looked at so far. 
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dN/d^ for WW production at ILC 
A=3 TeV, £=500 fb"' 




FIG. 14: The differential event rate dN/dcf), integrated over cos 6* as a function of for tlie 
SM and various mNCSM = 0) contributions with A = 3 TeV for W^W~ production at 
a 1 TeV ILC with 500 fb^^ integrated luminosity. Each set of NCSM data points is for 
c^j_i, = 1 for one choice of fiu, and all others being zero. The left column of plots is for 
K2 = 0, while the right column is for k,2 = 1- The top row of plots is for W^/^PVg~ 
production, while the bottom row is for Wf^W£ . The error bars, as before, include 
statistical error and luminosity uncertainty. 

2. dN /d cos e 

Instead of integrating over cos 6 and looking at dN/ dcj), we can also integrate over cf) and 
look at dN/ d cos 9 for different W helicity choices as a way of determining the parameters of 
the NCSM. In general, integrating over the full range of will remove any sensitivity to the 
NCSM contributions, since they typically vary sinusoidally with 0. Here we integrate from 
= to = TT and see that the cos 6 distributions in this half of the detector are sensitive 
to the NCSM parameters. 

For example. Fig. [T6l shows dN/dcos6 for W^W^ production for the half-cylinder defined 
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dN 



with cos6'g[-1,-0.9] for WrW^q production at the ILC 




FIG. 15: dN/d^ for W^Wq for the bin -1 <cos9 < -0.9. 

above. In this case, the NCSM distribution with cis and /t2 = 1 is clearly distinguishable from 
those of the SM and the other 11 NCSM parameter choices, all of which have negligible event 
rate for this helicity choice. As shown in Fig. [T71 on the other hand, for Wfl~W^ production 
the Ci2 = 1 and C13 = 1 with K2 = distributions show deviations from each other, as well 
as the SM. 



V. CONCLUSIONS 

The Noncommutative Standard Model is a model of physics with a rich history. Motiva- 
tion for its study comes both from theoretical considerations from string theory, and generic 
properties of quantum gravity, as well as from the NCSM's unique and rich phenomenology. 

In this paper we have contributed to the large body of work on the phenomenology of 
the NCSM by investigating WW production in this model. We have found that in both 
W+W- W+W- and e+e" W+W' in the NCSM, unlike in the SM, the cancellation 
of terms in the amplitude that grow with s and s^, respectively, does not occur. At some 
value of s, these terms violate partial-wave unitarity, and we have estimated that value. 

For an interesting range in NCSM parameter space (for high enough A and low enough s). 
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with (p<jT for W iW^h production at the ILC 
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FIG. 16: The differential event rate dN/d cos 9 integrated over the half-cyhnder from = 
to = vr for the SM and the NCSM {k,2 = 1) contributions for W^W£ production at a 
1 TeV ILC with 500 fb"'^ integrated luminosity. The data points represent the NCSM 
event rate for cis = 1 and all other c^y set to zero, and the dashed curve is the SM event 
rate. The error bars show the one-sigma statistical errors on the NCSM event rate. 

however, the model satisfies partial wave unitarity. We have seen that for much of this range, 
the model has observable signatures at LEP, the LHC, and the ILC. We have determined 
the region of parameter space excluded by observations of the WW production differential 
cross section at LEP, and found that it extends to much higher values of the NC scale A 
then are ruled out by partial-wave unitarity at LEP energies. For the LHC, we have found 
that the search reach obtained from deviations from the WW production differential cross 
section and invariant mass distribution does not extend beyond the region that is ruled out 
by unitarity. This shows that the LHC is not an optimal tool (at least using this observable) 
to probe this model. 

For the ILC, the search reach, using WW production, extends to A ~ tens of TeV, well 
beyond the unitarity bound on A at ILC energies. We have shown that the polarization 
capabilities of the ILC are essential in probing regions of parameter space (via the left-right 
asymmetry) that are unobservable in the unpolarized cross section because of cancellations 
in the amplitude. Furthermore, measurement of the final state W polarizations provides 
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with (p<n for W qW^r production at the ILC 

dcos6' 

K2=Q, A=2 TeV, and £^500 fb"' 
200 1 , , , 




-1.0 -0.5 0.0 0.5 1.0 

cos9 

FIG. 17: Like Fig. [161 except for VFg''iy^ and k,2 = 0. The data points with thick (thin) 
error bars are for C13 = 1 (ci2 = 1), with all other Cf^iy set to zero. The dashed curve is the 

SM event rate. 

sensitivity to certain parameters of the model, namely C12, C13, C23, and cqs, whose contri- 
butions to the cross section cancel in the sum over final state helicities. Assuming efficient 
determination of the W helicities is possible, we have given examples where these NCSM 
parameters at may be measured at the ILC for values of A up to several TeV. 

There is much more work that can be done in this direction. One could go beyond 
calculating differential cross sections and instead incorporate the amplitudes into an event 
generator. This would allow the effects of parton showering, detector simulation, and Stan- 
dard Model backgrounds to be taken into account. This has been undertaken for other 
NCSM processes at the LHC by Alboteanu et al. j52] 

APPENDIX A: AMPLITUDE SQUARED FOR e+e'^W+W- WITH K2 = 

In this and the next appendix we present the helicity amplitudes for e~^e~—>-W~^W~ pro- 
duction in the NCSM with ^2=0. The next section we will present the K2 coefficients in 
the contribution to the helicity amplitudes from the NCSM with K2 7^ 0. The amplitude 
for e~^e~ —>■ W~^W~ in the NCSM has contributions from the Feynman diagrams given in 
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Eqs. nil. 121 IIII.13t and IIII.141 so that (for down-type fermions) 

Me+e-^VF+VK- = M'^^ + M|f + Ml"" + M^^^ + M^^ + Mj^^ + M,^^, (A.l) 

where M^^, M^^, M/^, and refer to the s-channel photon exchange, s-channel Z boson 
exchange, t-channel neutrino exchange (for down-type fermions this would be u-channel neu- 
trino exchange) and contact interaction diagrams, respectively. M^^ refers to the standard 
model contribution of a diagram, whereas M^*-^ refers to the 0(6') NCSM contribution. 
To calculate the cross section, we must square the amplitude. Doing this yields, to 0(^), 



where for a pair of diagrams a, 6, M is defined as 

M, 6 = 2Re(Mf ^Mfc^^ + M^^M^^) . (A.3) 

In this appendix we give the values of M for every diagram pair and every helicity choice, 
for K2 = 0. The helicities are specified in the order M{he-he+hw+hw-)- 

It is convenient to group the expressions by diagram pair and by initial state helicities. 
So, for example. Section A.l gives the values of M^,^, = 2Re(M^^^M^^ + M^^^M^^) with 
K2 = for a left-handed electron and a right-handed positron-the values are then listed by 
the helicities of the final-state W bosons. The prefactor that all the expressions share is 
given first, while the remaining factor for each W helicity choice follows. So the first line in 
Section A.l, 

00= se{-c^9o2 + s^9oi){3 + f3^-4p^)s, (A.4) 

gives the value M^^-y^ (LROO) if multiphed by the prefactor given at the beginning of Section 
A.l. If any polarization choice is not shown, it is because the corresponding contribution is 
zero. 

1. Values of M^^^^{LRhy^+h-i^-), with K2 = Q 

Prefactor: g^Q'fS^f3s/ (8m^) 
00= se{-c^eo2 + s^eoi){3 + (3^ - 4f3^)s 

OL = 4mlr{{-2 + sl + 2ce)e,2(3 + (1 - ce)ses^(^oi - dnP) + (-1 + ce)s,c^(^02 - ^23/?)) 
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OR = 4m^((-2 + sl- 2ce)9i2p + (-1 - ce)sec^{9o2 + ^23/9) + (1 + ce)sgs^{9oi + 
LO = 4m^((2 + 2ce - s^Oup + (1 + ce)sec^(-^02 + ^23/9) + (1 + ce)ses^(^oi - ^13/?)) 
LL = -4m^((-2 + + 2ce)^i2/9 + (1 - Q)sec^(^o2 + ^23/?) + (-1 + ce)ses<^(^oi + ^i3/9)) 



2. Values of My^-y^(i2L/i^;j/+/i^;j/-), with K2 = 

Prefactor: g'^Qjs^wPs/ {Sm^,) 
00= -se(-cA2 + Voi)(3 + /5^- 4/5=^)5 

OL = 4m^((2 + 2cg - sl)9uP + (-1 - ce)s9S^(^oi + ^i3/9) + (1 + ce)sgc^{9o2 + ^23/9)) 
Oi? = -4m^^.((~2 + si + 2ce)^^i2/? + (1 - cg)ses^{9oi - 9nf3) + (-1 + ce)sec^{9o2 - ^23/^)) 
LO = 4m^((-2 + sl + 2ce)9i2l3 + (1 - ce)sec^(^02 + ^23/9) + (-1 + ce)sgs^{9oi + ^i3/9)) 
LL = 4m^((-2 + - 2c0)^i2/5 + (1 + ce)s0C^(^o2 - ^23/5) + (1 + C0)seS0(-^oi + ^13/?)) 



3. Values of Mz^j^iLRhw+hw-), with K2 = 

Prefactor: (-c/;^ - gl,)g*QfS^^(3s/ (32c^m^(s - m|)) 
00 = (-1 + 2slr)se{-c^9Q2 + sAi)(3 + - 4/5=^)5' 
OL = -2m^((l - c,)s,c^(^02(-3 + 4,s^,. + /J^) + 2 (l - 25^^) ^^23/9) 

+ (-1 + ce)ses^{9^^{-Z + 45^^ + + 2 (l - 2s^) ^13/?) 

- 2 (2s^ - 1) (-2 + + 2ce)^i2/?)s 

OR = 2m^((-l - ce)sec^(-3^02 + 4^(^02 + ^23/9) + ^02/^' - 2^^23/9) 

+ (l + c,)s,s<^(-3^oi + 4s2^(^oi + ^i3/9) + ^oi/9'-2^i3/3)-2(2s2^-l) (2 + 2c,-s2)^i2/3)5 

LO = 2m^((-l - ce)sec^(^02(-3 + 4s^ + z^^) + 2 (l - 25^^) ^23/9) 

+ {l + ce)s6S^{9oi{~?, + Asl, + l3^) + 2 (l-2s^) ^i3/3) + 2 (2s2^-l) (2 + 2c0 -s^)^i2/3)s 

LL = 8sem^(c<^6'o2 - s<^6'oi - c^92zl3 + s^9i2,l3) 

LR = -2m^((l - ce)sec^(-3^02 + 45^(^02 + ^23/?) + ^02/9' - 2^23/9) 

+ (-1 + ce)sgs^{-?,9oi + 45^(^01 + ^13/9) + ^01/9' - 2^13/9) 

+ 2 (2s^ - 1) (-2 + + 2ce)^i2/9)s 

i?0 = 8semt^(c^(^02 + ^23/9) - 5^(^01 + ^13/9)) 



4. Values of M^^^^ (i?L/ivF+^Ty-)? with K2 = 

Prefactor: - gi)g'^Q js^^ps/ (32c^m^(s - m|)) 
00 = (-1 + 2sl,)sg{-c^9o2 + s^eoi)(3 + /9^ - 4/3')? 
OL = 2m^((-l - ce)sec^(-3^02 + 4s^(^o2 + ^23/9) + ^02/9' - 2^23/?) 

+ (l+ce)seS0(-3^oi+4s^(^oi+^i3/9)+^oi/9'-2^i3/3)-2(2s^-l) (2+2ce-s^)^i2/3)s 
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OR 
LO 

LL 
LR 

RO 



-2m^((l - ce)s0C^{0o2{-S + + Z^^) + 2 (l - 2s^) 023^) 

+ (-1 + ce)ses^{eoi{-3 + Asl, + p^)+2{l- 2s^) ^13^) 
- 2 (2s^ - 1) (-2 + sl + 2ce)e,2P)s 
-2m^((l - ce)s0C^{-3eo2 + 4s^(^^o2 + ^^23/9) + ^^02/?' - 2623^) 

+ (-1 + ce)s,s^(-3^oi + 4s2^(^^oi + ^^13/3) + ^^01/9' - 2^^13/9) 

+ 2 (2s^ - 1) (-2 + + 2ce)euP)s 

8sem^(c^(^02 + ^23/3) - s<^(^oi + ^is/^)) 
2m^((-l - ce)sec^ieo2{-3 + + ^2^, ^ 2 (l - 25^^) ^23/9) 
+ (l + ce)ses<^(^^oi(-3 + 4s^ + /52)^2 (l-2s^) ^i3/3) + 2 (2s^-l) (2 + 2ce -s^)^i2/3)s 
8s0m^(c</,6'o2 - S06'oi - c^623(^ + s^O^P) 



5. Values of Mz^ZsiLRhw+hw-), with K2 = 



Prefactor: -(51;^ + gyfg'^Ps'^/ {QAc^m'^is - m\f) 
00 = " 

0L = 
0i? = 



LO 

LL 
LR 



2^w'((~l + C0)s0C^(6'o2(-l + 2s 



+ (-1 + C0)ses^{9^ 



01 



2s 



+ P^)-2s 

- O13/3) 



W^23P) 
2\ 



2slri-2 + sl + 2ce)ei2l3)s 

1 - ce)sec^{-eo2 + 2s^(^o2 + ^23/9) + ^02/9l 
+ (1 + ce)s,s<^(-eoi + 2s^(^oi + ^13/5) + ^oi/5') + 2s^(-2 + s^ - 2ce)eu(3)s 
2m^((l + C0)sec^(^O2 - 25^(^02 - O23P) - ^02/3') 

+ (1 + ce)ses^(eoi{-l + 2s^ + p^) - 25^^13/?) + 2s^(2 + 2ce - s^)^i2/9)s 
8se"7,^(c^^02 - s^9oi - 0^623^ + S^9i3(3) 
-2m^((l - ce)sec^(-^02 + 2s2^,.(^o2 + ^23/9) + ^02/9') 
+ (-1 + c,)s,s^(-^oi + 2s^(^oi + ^13/9) + ^01/?') + 2s^(-2 + + 2ce)^i2/?)s 



i?0 = 8S0m|j.(c^(^o2 + ^23/9) - 5^(^01 + ^13/9)) 



6. Values of Mz^Zs{RLhy^+hy^-), with K2 = 



Prefactor: -{g{ - glfg^(5s^l {Uc^ml,{s m\f) 



00 
OL 

OR 

LO 

LL 
LR 



-sivSg{-C^e02 + sAl)(3 + P~ 
-2m^((-l - Ce)s,C^(-^02 + 2S^(^02 + ^23/9) + ^02/9' 



+ (H-ce)ses^(-^oi + 2s^( 



-2m^^((-l + c,)sec^(^^02(-l + 2s^^ + Z^^) - 2s2^^^23/3) 



?01+^13/3)+^Oir 

,2 



+ (-1 +C0)ses^(6'oi 



2s2vk(^^oi - h3^) - OoiP') 
2m^((l - ce)sec^i-eo2 + 2s'^(^^o2 + ^23/9) + ^02/9': 
+ (-1 + ce)ses^{-eoi + 2s2,(^oi + O13P) + ^oi/9') 
8sem^(-c<^(^02 + O23P) + s^iOoi + 9i3P)) 
2m^((l + cg)sec^(eo2(-l + 2s%^ + 0") - 2s^^23/9) 

,2 (a a a\ a /32\ , r,„2 



2s^(- 



2s^( 



2 + s^ - 2cg)ei2P)s 
2 + sl + 2ce)ei2P)s 
-2 + sl + 2ce)ei2p)s 



+ (1 + ce)ses<^(eoi - 25^(^01 - - Ooin + 2s^k(-2 + - 2cg)euf3)s 
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7. Values of Mf^y_^{LRhiY+h\y-), with K2 = 

Prefactor: g^QfS^^rs/ (64tm^) 
00 = -se{-c^9o2 + s^9oi){-l+P^){P^ + 2ce - 3/3)s' 
OL = _s(-2m^(-ei2/3(-2(/3 + 1)^ + 2ce((/3 + 1)' - s^) + s^S + /3(2 + 

+ sec^{2sl -{(3 + 1)2 + ce(/5 + l)2)(-^02 + ^23/9) 
+ ses42sl -{(3 + ly + ce(/5 + iy){doi - ^^la/?)) 
+ se(-c^(^02 + ^23) + s^(^^oi + 0is))f3{-l +ce + i-l + ce)P^ - 2(c, - l)ce/?)s) 
Oi? = s{-2ml,{euPW - If + 2c,((/5 - If - si) - sl{3 + (-2 + f3)P)) 

- ses4i(3 - If - 2sl + ce{P - iffOoi + 6^^^) 
+ sgc^{{P - If 2s] + ce{P lf ){e^^ + 92,(3)) 

+ Se(-C<^(^02 + ^23) + 50(^01 + ^13))/3(1 + Cg{P - If + /3(-2 + /3 + 252))5) 

LO = s(2rr4(^i2/3(2(/3 - 1)^ + 2c,((/3 - 1)^ - s]) - s^(3 + (-2 + 

+ s,c^((/3 - If - 2sl + ce{(3 - lf){-9o2 + ^23/9) 
+ ses^{{(3 - If - 2sl + ce{(3 - lf){9^^ - 9^^(3)) 
+ seM92z - 9o2) + s^{9oi - 9n))P{l + ce{(3 - if + (3{-2 + /3 + 2sl))s) 

LL = 4se{-c^{9o2 + 05^23) + ^^(^oi + C0^i3))m^/3(1 + - 2ce(3)s 
LR = 16se{-2 + - 2cq){-c^9q2 + 5^6*01)^^4' 

RQ - s{2ml,{9,2(3{-2{(3 + 1)^ + 2c,((/3 + 1)^ - s^) + s^(3 + /?(2 + /?))) 

- sec^(2s2 - (/5 + 1)2 + ce(/5 + 1)2)(^02 + 6^23/9) 
+ ses^(2s2 - (/3 + 1)2 + Ce(/3 + l)2)(^^oi + 6^13/9)) 

- se(c^(^23 - ^02) + s^{9oi - ^i3))/9(-l + ce + (-1 + ce)/?^ - 2(ce - l)ce/3)s) 
i?L = -16s0(-2 + sl + 2ce){-c^9o2 + s^9oi)m^ 

RR = -4se(c</,^02 - 5,^^01 - C0C^923 + ces ^9 13) 771^(3(1 + - 2ce(3)s 

8. Values of Mf^Zs{LRh\Y+h\Y-), with K2 = 

Prefactor: {-g{ - g(.)g'^s^/ (128c^tm^(s - m|)) 

00 = -slySe{-c^9o2 + s^9oi)i-l + p^)ip'' + 2ce - 3p)s^ 

OL = -s(m2^(-ses^(2s2 - (/3 + 1)^ + ce(/3 + l)2)(6^oi(-l + 2sl, + Z^^) - 2slr9n(3) 

+ sec^{2sl -{(3 + If + ce{(3 + 1)2)(^02(-1 + 2sl, + Z^^) - 2s^^23/5) 
+ 2s^^i2/5(-2(/3 + 1)2 + 2c,((/3 + 1)2 - s2) + ^2(3 + /3(2 + /?)))) 

- C^Se(-C0(^O2 + ^23) + S^{9oi + 9i3))P{-l +C0 + (-1 + Cg)P^ - 2{cg - I)c0f3)s) 

OR = s(m^(-s,c^((/3 - 1)2 - 2^2 + c,(/3 - l)2)(-^02 + 2^2,(^02 + ^23/?) + ^02/?') 

+ ses^{{(3 - 1)2 - 2^2 + ce{(3 - l)2)(-^oi + 2^(^01 + ^13/?) + ^oi/5') 
+ 2s2,^i2/9(-2(/? - 1)2 + 2c,(s2 -{13- 1)2) + s2(3 + (_2 + /5)/5))) 

- cl^S0{-c^{9o2 + 923) + s<^(^oi + 9n))P{l + ce(/9 - if + P{-2 + P + 2sl))s) 
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LO = s{ml,{-S0C4,{{P - If - 2sl + Ce{l3 - 1)')(^02(-1 + + 13^) - 2s^^23/?) 

+ ses^{{(3 - If - 2sl + ce{(3 - l)')(^oi(-l + 2s'^ + (5^) - 2s^^i3/3) 
+ 2s^^i2/5(2(/5 - 1)' + 2ce{{fi if s]) - 5^(3 + (-2 + (3)m 
- cl,se{c4,{92^ - 9o2) + 5^(^01 - ^i3))/3(l + ce{(3 - if + /3(-2 + /3 + 2s^))s) 
LL = 4s5im^(ce(m^(2c<^6'o2 - 2s<^6'oi - 20^6*23/3 + 2s<^6'i3/3) 

+ C^/?(C^(^23 - 2^02/5 + ^23/9') - S^(^^13 " 2^^01/3 + ^^13/3'))s) 

+ /5(2m^(-c^^o2 + s^(^oi - ^13/3) + c^^23/5) 

+ C^(C<^(^02 + ^^02/?' - 2c2^^23/3) - 5,^(^01 + ^Ol/S' " 2c2^i3/3) )s) ) 

Li? = 8 (2s^ - 1) se(-2 + - 2ce) (-c^^02 + S0^oi)m^ 

m = s(m^(-sec^(2s^ - (/3 + 1)' + ce{(5 + 1)') (-^02 + 2s^(^o2 + ^23/?) + ^02/?') 

+ ses^(2s^ - (/3 + 1)2 + ceiP + l)')(-^oi + 2s^(^oi + ^i3/3) + ^oi/3') 
+ 2s^^i2/3(-2(/3 + 1)' + 2c,((/3 + 1)2 - s^) + s^(3 + /3(2 + /3)))) 
+ c^se(c^(^23 - ^02) + 5^(^01 - 0i^))(3{-l + ce + (-1 + c^)/?' - 2(c0 - l)ce/3)s) 
RL = -8 (2s^ - 1) se{-2 + + 2ce) (-c^^02 + s^O^Ww 
RR = -As0m^{c0{2m^{-c^{eo2 + ^23/9) + s</,(^oi + Onf^)) 

+ (3{2m^w{c^{eo2 + 023^) - s^iOoi + OisP)) 

+ C^(-C^(^02 + ^02/3' + 24623(3) + S^idoi + ^01/?' + 2c^^i3/?))5)) 



9. Values of Mf^y^{LRhy^+hy^-), with K2 = 

Prefactor: g^QfS^s/ (64-um^) 
00 = s,(-c^^^02 + +/^')(/^' - 2ce - 

OL = s{2mlr{e,2(3{2{(3 - if + 2ce{sl -{(3- if) - s^(3 + (-2 + (3)(3)) 

+ sec^{2sl - (/? - 1)^ + ce{P - l)')(-^02 + ^23/9) 
+ s,s^(2s2 - (/j - 1)2 + ce{(3 - lf){eoi - e,s(3)) 
+ se{c^{e23 - 602) + s^iOoi - e,,))(3{l - ce{(3 - if + (3{-2 + /? + 25^))^) 
OR = -s(2m^(^i2/3(2(/3 + If + 2c0{(P + if - sj) - s^(3 + /3(2 + P))) 

- ses^{{(3 + If - 2sl + C0{(3 + l)2)(eoi + ^13/9) 
+ sec^{{(3 + If - 2sl + ce{(3 + 1)')(^02 + ^23/3)) 
- s,(c^(^23 - ^02) + s^{eoi ~ ^i3))/5(l + c,(/3 + If + p(2 + P- 2sl))s) 
LO = s(2m^(^i2/?(2(/3 + 1)^ + 2ce{{(3 + if - s^) - s^(3 + (3{2 + (3))) 

+ sec^iiP + 1)' - 2sl + ce{(3 + lf){-Oo2 + 023(3) 
+ ses^((/3 + 1)' - 2sl + ce{(3 + l)2)(^oi - ^^i3/9)) 
+ se(-c^(^02 + ^23) + s^{eoi + ^i3))/9(l + c,(/3 + 1)2 + (3{2 + (3- 2sl))s) 
LL = 4s,m^/3(s^(^oi(l + + 2c,/3) - ^13(0^ + 2/? + c,/?^ - 2s^/5)) - 0/02(1 + + 2ce/3) 

+ c^^23(ce + 2/3 + - 2^2/3) )s 

LR = 16s0(-2 + - 2ce)(-c/o2 + s^^'oO^^vk 
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RO = s(2m^(ei2/3(-2(/3 - 1)^ + 2c,((/3 - if - + s^(3 + (-2 + /?)/?)) 

- s,c^(2s^ - (/3 - If + cei/3 - \f){Q^2 + ^23/9) 

+ se(-c^(^02 + ^23) + 5^(^01 + ^i3))/9(l - ce(/? - 1)' + /3(-2 + /3 + 2s^))s) 
i?L = -16se(-2 + + 2ce)(-c<^6'o2 + s<^6'oi)"^vk 
RR = 4sem^/5(-c^(^o2(l + + 2c,/5) + 623(09 + 2(3 + ce(3^ - 2slf3)) 

+ s</,(^oi(l + 13^ + 2cel3) + Oi^ice + 2(3 + eels'" - 2sl(3)))s 



10. Values of Mf^ZsiLRhw+hw-), with K2 = 

Prefactor: {—g^ — 9v)g'^''^'^/ (128c^Mm^(s — m%)) 
00 = sl,se{-c^9o2 + s^eoi){~l + P^){P' - 2c, - 3P)s^ 

OL = s{ml,{-sec42sl (/? - 1)^ + ce(/? - 1)')(^^02(-1 + 25^^ + /3^) - 2s2p,.^^23/3) 

+ ses^{2sl -{(3- If + ce{P - l)')(^oi(-l + 2^ + - 25^^13/?) 
- 2s^ei2/3(-2(/5 - If + 2c,((/5 - 1)^ - s]) + sl{3 + (-2 + 

- cl,seMe2s - 602) + s^iOoi - e^s))f3{l - ce{(3 - if + (3{-2 + (3 + 2sl))s) 
OR = -s(mlr(-seS4,((P + if - 2sl + ce{(3 + l)')(-0oi + 2s^(^oi + ^i3/9) + ^oi/5') 

+ sec^iiP + If - 2sl + ce{f3 + l)')(-^02 + 24.(^^02 + ^23/?) + ^02/?') 
+ 2s^^i2/?(2(/? + If + 2ce{{f3 + if - s]) - sl{2> + (3(2 + (3)))) 
+ c^se(c^(^23 - ^02) + 5^(^01 - 0^3))P{1 + ceiP + If + P{2 + P- 2sl))s) 
LO = s{ml,{-sec^{{f3 + if - 2sl + ce{P + 1)')(^02(-1 + 2s^ + (3^) - 2sl,e2-il3) 

+ seS4,{{P + If - 2sl + ceiP + l)')(^oi(-l + 2s^ + (3^) - 2s^^^i3/5) 
+ 2sl,e,2mP + If + 2ce{{P + If - si) - 5^(3 + /3(2 + 

- clrSe{-c^{eo2 + 023) + ^^(^01 + ^i3))/9(l + ce{P + If + (3{2 + (3- 2sl))s) 
LL = As0m'^{c0{2rrf^{c^9o2 - s^Oqi - €^623^ + s^9i3(3) 

+ C^l3{-C^{e23 - 2^02/3 + ^23/3') + 50(^13 - 2^01/3 + ^13/3'))S) 

+ P{mlr{2c^eo2 - 28^001 - 20^623(3 + 2s^ei3P) 
+ c^(c^(^02 + ^02/3' - 2c^^23/3) - s^ieoi + ^01/9' - 2clenP))s)) 
LR^ 8 (2s^ - 1) se(-2 + - 2ce)(-c^eo2 + Voi)^^ 

RO = s{m^^{-sec^{2sl -{13- if + ce{P - l)')(-^02 + 25^,(^02 + ^23/9) + ^02/9') 

+ ses^i2sl -{13- If + ce{(3 - l)')(-^oi + 2sl,{Q^^ + ^13/3) + ^oi/3') 
+ 2s^^i2/3(-2(/3 - If + 2ce((/3 - 1)^ - si) + sl[3 + (-2 + /3)/3))) 
- c%,se{-c^{B^2 + ^^23) + 5^(^01 + ^i3))/3(l - ce{{3 - if + (3{-2 + (3 + 2sl))s) 
RL = -8 (2s^ - 1) se{-2 + + 2ce)(-c^^o2 + s^^w 
RR = 4sem^(ce(2m^(c^(^o2 + ^23/3) - s^(^oi + ^i3/3)) 

+ C^/3(C^(^23 + 2^02/3 + 023(3^) - S^{e^3 + 2^01/3 + ^13/3'))S) 

+ /3(2m^(c<^(^02 + ^^23/3) - 5,^(^01 + ^i3/3)) 
+ 4(c^(^o2 + ^02/3' + 2c2^23/3) - 5,^(^01 + ^01/3' + 2cle,3(3))s)) 
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11. Values of Mf^f^{LRh-^+h-^-), with ^2 = 



Prefactor: g'^sgs^/ (512t^m^) 
00= (--c<^(^o2 + ^23)-S0(^oi + ^i3))((/3+l)'-ce(l + /3(-4s^ + (2 + /5)(2 + /3(2 + /3)))) 

~ 2sl{l + 13 + (3^ + 3p^))s 

OL = (-c^(^02 + 023) + 5^(^01 + ^i3))((/? - 1)' + ce(l + /?(4s^ + (-2 + /3)(2 + (-2 + (3)f3))) 

+ 2sl{-l + /3 + /3^ -3/3^))s 

OR = (c^(^23 - ^02) + s^{9oi - 9n)){{P - + cg{l + P{4sl + (-2 + (3){2 + (-2 + (3)(3))) 

LO = 8(-c^(^o2 + ce^23) + s^iOoi + ce^i3))m^(ce + SceP'' - /3(3 + - 2s^)) 
LL = -8(-c^^o2 + 5^^oi)m^(2(/? - 1)^ + 2ce(l + f3{-2 + + (3)) - sl{l + (-4 + f3)f3)) 
LR = (-c^(^23 - ^02) - 8^(601 - eisMP + 1)' - ce{l + P{-^sl + {2 + (3){2 + (3{2 + (3)))) 

-2sl{l+(3 + (3^ + ?>(3^))s 
RO = 8(-c/o2 + s^9oiWw{2{P + l)^ - 2ce{l + {3{2 + P- sD) - sl{l + /5(4 + /?))) 
RL = -8{c^9o2 - s^9oi - C0C^923 + Cgs^9i3)m^{ce + Sc^/?^ - (3{2> + (3^ - 2s^)) 



12. Values of Mf^f^{LRhw+hw-), with K2 = 

Prefactor: g^sgs'^/ (512-u^m^) 
00 = (c^(^23 - 9o2) + s^{9oi - 9i3)){-{P - 1)^ + ce{l + /3(4s^ + (-2 + P){2 + (-2 + P)/3))) 

- 2sl{-l + f3 + (3^ - 3(3^))s 

QL = {-c^{923 - 9o2) - s^{9oi - ^i3))(-(/i^ + 1)' - c,(l + (3{-Asl + (2 + f3){2 + /i(2 + /?)))) 

+ 2s^(l + /3 + /33 + 3/32))s 

Oi?= ( c^(^02 + ^23) -5^(^01 +M)(-(/?+l)'-ce(l + /5(-4s^+(2 + /5)(2 + /?(2 + /?)))) 

+ 2s^(l + /3 + /3^ + 3/32))s 
LO = -8(c<^^02 - 5,^^01 - cec^923 + ces^9i3)rn^{c0 + 3ce/3^ + /3(3 + - 2s^)) 
LL = 8(-c^eo2 + s^^oi)m^(-2(/3 + 1)^ - 2ce(l + /3(2 + /3 - s^)) + + /3(4 + 
Li?= (_c^(^o2 + ^23) + s^(^oi + ^i3))(-(/3-l)' + ce(l + /3(4s^ + (-2 + /3)(2 + (-2 + /3)/3))) 

-2s2(-l+/? + /?3_3/52))s 

/?0 = -8(-c<^^02 + s^9m)ml,{-2{P - if + 2ce(l + /3(-2 + + /?)) + + (-4 + 
= 8(-c^(^o2 + ce^23) + s<^(^oi + ce9iz)Ww{ce + 3ce/3' + (3{2> + (3^- 2sl)) 
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13. Values of Mcy^{LRh^r+h■^Y-), with K2 = 



Prefactor: g^Qfs'^sgPs/ (IGmly) 

00 = (1 - ce){c^62z - s^6i3 - c^6q20^ + s^6qiP'^)s 

OL = (1 + Cg){c^623 - S<p6i3 - 0^602^"^ + s^6qiP'^)s 
0R= (1 + ce){-c4923 + eo2P^) + s^iOn + 9oi(5^))s 
L0= -4:{ce-l){-c^e23 + s^ei3)ml,(3 

LL^ (1 - C,)(-C^(^23 + 002^^) + S^{ei3 + eoiP^))s 

LR = 4(c0 - 1) (-C0^23 + s^ei3)w?w(5 

14. Values of McZs(-C'-R^vf+^vf-)! with K2 = 
Prefactor: {g{ + g(,)g^S0(3s'^ / (32m^(s - m|)) 

00 = (1 - Cg){c^923 - S^0i3 - C^6o2f3'^ + s^6oi(3'^)s 

OL = (1 + ce){c^623 - s^9i3 - c^9q20^ + s<^^oi/3^)s 
Oi? = (1 + ce) (-c^(^23 + ^02/3') + s^{ei3 + ^oi/5'))s 
LO = -4(ce - 1) (-C0^23 + s^9i3)m^w(5 

LL^ (1 - C) (-0^(^23 + ^02/?') + 5^(^13 + ^01/9'))S 

LR = 4(c9 - l)(-c^^23 + s^ei3)m'^P 



15. Values of Mcft{LRhw+h\Y-), with At2 = 

Prefactor: ^S^sV (128tm^) 

00 = (2s^ - (/3 + 1)2 + c,(/3 + l)')(c^e23 - s^e,3 - C^^02/3' + s^9oiP'')s 

OL = - 1)2 - + ce(/5 - l)2)(c^^23 - s^e,3 - c^^02/9' + s^eoi(3^)s 

OR = {{(3 - 1)2 - 2s2 + ceilS - l)2)(-c^(^23 + ^02/?') + s^(^i3 + ^oi/?'))s 

LO = -8(-c^^23 + s^ei3)mly{ce - 

LL = (2^2 - (/? + 1)2 + C0{(3 + 1)2)(-C^(^23 + ^02/3') + 5^(^13 + ^01/9'))s 

Li? = 8(-c^^23 + s^ei3)m^w{ce - 



16. Values of Mcf„{LRh^^+hy^r-), with K2 = 

Prefactor: g^sgs"^/ (128-umL) 

00 = (2^2 -{p- if + ceiP- l)')(c/23 - s^e,3 - C0^o2/9' + s^eoiP^)s 

OL = {{(3 + 1)2 - 2s2 + Cg{l3 + 1)2)(C^023 " S<^^13 " 0^2/^' + S0^^Ol/9')s 

Oi? = {{/3 + 1)2 - 2s2 + ce{/3 + l)2)(-c^(^23 + ^02/?') + 5^(^13 + ^oi/3'))s 

LO = 8(-c/23 + s^e^3)mlr{(3^ - cl) 

LL = (2s2 -(13- 1)2 + ce{[3 - l)=^)(-c^(^23 + O02P') + ^0(^13 + OoiP'))s 

LR = 8(-c^^23 + S0^i3)m^(c^ - /?2) 



37 



APPENDIX B: AMPLITUDE SQUARED FOR e+e'^W+W-, K2 TERMS 



In this appendix, we give the values of the coefficients of K2 in the NCSM contributions 
to M for every diagram pair and every hehcity choice. The notation and organization is the 
same as in the previous appendix. For example, section B.l gives the @{k2) contribution 
to the quantity M^^^^ = 2Re(M^J^M^^ + M^/^M^^) for a left-handed electron and a 
right-handed positron. The line 

00= -se(-cA2 + sAi)(-3 + /9')(/9'-l)'s' (B.l) 

gives the coefficient of k,2 in My^^^ (Li?00) , if this expression is multiphed by the prefactor 
given at the beginning of Section B.l. If any polarization choice is not shown, it is because 
the corresponding contribution is zero. 



1. Coefficients of K2 for M~^^j^{LRhy[r+hy^-) 

Prefactor: g^^Qjs^wP^/ (Sm^^) 

00 = -sg{-c^9o2 + s/oi)(-3 + - if 

OL = -Amlr{-1 + /3')((-2 + sl + 2ce)9uP + (1 - ce)sec^{eo2 + 023/3) 

+ (-1 + ce)ses^(^oi + ^i3/9))s 

OR = -4m^(-l + /3')((-2 + - 2ce)^i2/3 + (1 + ce)secjeo2 - ^23/9) 

+ {l + ce)ses^{-eoi + ei3l3))s 

LO = 4m^(-l + /3')((-2 + - 209)612^ + (-1 - ce)sec^{eo2 + ^23/?) 

+ {1 + ce)ses^{0oi + enP))s 

LL = 32s0m^{c^9o2 - s^Oqi - c^92z(3 s^9i2,(3) 

LR = 4m^(-l + /3')((-2 + + 20^)^12/3 + (1 - ce)ses^(^oi - ^i3/9) 

+ (-1 + C0)s0C^{9o2 - 92zP))s 

RO = 32s0m^(c^(^o2 + ^23/?) - s^(^oi + ^13/9)) 



2. Coefficients of K2 for M~^^j^{RLhy[r+hy^-) 

Prefactor: g^^Qjs^wP^/ (Sm^^) 

00 = S0{-c^9o2 + Voi)(-3 + P^) {P^ - if 

OL = 4m^(-l + /3')((-2 + - 2ce)^i2/5 + (1 + ce)sec^{9o2 - ^23/?) 

+ (l + ce)ses^(-^oi + ^i3/9))s 
OR = 4m2j,(-l + /32)((-2 + sl + 2cg)9uP + (1 - ce)s,c<^(^02 + ^23/9) 

+ (-1 + C0)S0S^(^O1 + ^13/3))S 
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LO = -4m^(-l + /3')((-2 + si + 2ce)9i2P + (1 - ce)ses^{9oi - OisP) 

+ (-1 + C0)s0C^{6o2 - 023P))S 

LL = 32s9m^(-c<^(^02 + ^23/?) + s<^(^oi + OisP)) 

LR = -4m^(-l + /3')((-2 + - 2q)^i2/3 + (-1 - ce)sec^{eo2 + 02^(3) 

+ (1 + ce)ses^{eoi + ^i3/?))s 

i?0 = 'i2serr4v{-c^6o2 + s</,6'oi + C06'23/5 - s^6is(3) 
3. Coefficients of K2 for Mz^-yX^Rh-^+h-^-) 

Prefactor: ((?^ + gi)g'^Qfsl. (2s^ - Ss^^ + 1) /3sV (32c^m^(s - m|)) 
00 = -s,(-c/o2 + s^9oi){-3 + - 1)' 

OL = -4m^(-l + P^){{-2 + si + 2ce)ei2p + (1 - ce)sec</,(^02 + ^23/5) 

+ (-l + ce)ses^(^oi + ^i3/9))s 

Oi?= 4m^(-l+/32)((2 + 2c0-s^)^i2/3+(l+ce)sec.^(-^o2 + ^23/9) + (l+ce)s,s<^(^oi-^i3/3))s 
LO = 4m^(-l + /3')((-2 + - 2^)^12/? + (-1 - c,)s,c^(^o2 + ^23/9) 

+ {1 + ce)s0S4,{eoi + ei3P))s 

LL = 32sem^(c06lo2 - 5^6101 - c^023/3 + s^9i2,P) 

LR = 4m^(-l + /3')((-2 + + 2ce)^i2/3 + (1 - ce)s0S</,(^oi - Or,P) 

+ (-1 + ce)sec^{9o2 - ^23/3))s 

i?0 = 32s0m^(c^(^o2 + ^23/3) - 5^(^01 + ^i3/9)) 



4. Coefficients of K2 for Mz^-yX^^^w+^w-) 

Prefactor: (fir^ - g{,)g^Qfsl, (2sfy - 3^2, + 1) (32cVmt^(s - m|)) 
00 = -se(-c<^^02 + s^^oi)(-3 + - if 

0L= 4m^(-l+/?2)((2 + 2ce-s^)^i2/3+(l+ce)sec^(-^02 + ^23/9) + (l+ce)ses<^(^oi-M)s 
Oi? = -4m^(-l + /3')((-2 + + 2ce)^i2/9 + (1 - ce)sec^(^02 + ^23/9) 

+ (-1 + C0)sgS^{eoi + ^13/9))S 

LO = 4m^(-l + /3')((-2 + + 2c,)^i2/3 + (1 - c,)s,s^(^oi - OnP) 

+ (-1 + C0)sgC^{9o2 - 923(3))S 

LL = 32sem^(c^(^o2 + ^23/?) - s^(^oi + ^i3/9)) 

Li? = 4m^(-l + /?')((-2 + - 2ce)^i2/9 + (-1 - ce)sec^{eo2 + ^23/?) 

+ (l + ce)ses^(^oi + ^i3/5))s 

i?0 = 32sem^(c06'o2 - s^9qi - c^923/3 + s^9i2,(3) 
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5. Coefficients of K2 for Mz^ZsiLRh^^+h^r-) 



Prefactor: {g{ + g{,fg^sl,(5s^/ {Ucl,m%,{s - m\f) 
00 = se{-c^e^2 + s<^^oi)(-3 + f) {0" - if 

OL = 4m^(-l + /?')((-2 + + 2c0)euP + (1 - ce)sec«i(^^02 + ^23/9) 

+ (-1 + ce)ses^{eoi + ^i3/3))s 

OR = 4mlr{-l + /?')((-2 + sl- 2ce)eu(^ + (1 + c,)s,c<^(^02 - ^23/9) 

+ (1 + ce)s9S^(-^oi + ^i3/9))s 

LO = -4m2^(-l + /9')((-2 + - 2ce)^i2/3 + (-1 - c,)s,c^(^o2 + ^23/9) 

+ (1 + ce)s0S^(eoi + eiMs 

LL = 32s0m^{-c^9o2 + s^Oqi + c^92z(3 - s^OisP) 

LR = -4m^(-l + /5')((-2 + + 2ce)^i2/5 + (1 - ce)ses^(^oi - ^13/?) 

+ (-1 + C0)s0C^{6o2 - 6'23/9))s 

i?0 = 32sem^(-c^(^o2 + ^23/9) + s^iOoi + ^13/?)) 



6. Coefficients of K2 for Mz^ZsiRLih^r+hiY-) 

Prefactor: {g{ - glfg^s-^^^s^ / {Ucl,m'^{s - mlf) 
00 = -S0{-c^e^2 + s^^oi)(-3 + - if 

0L= 4m2^(-l+/52)((2 + 2ce-s2)^i2/9+(l+c,)s,c<^(-^02+^23/9) + (l + ce)ses<^(^oi-^i3/5))s 
Oi? = -4m2^.(-l + /?2)((-2 + + 2c0)9i2p + (1 - ce)sec^(^^02 + ^^23/^) 

+ i-l + C0)s0S^{eoi + ei3p))s 
LO = 4m^(-l + /3')((-2 + + 2c,)^i2/3 + (1 - 00)808^(601 - ^^i3/5) 

+ (-1 + C0)S0C^{6q2 - 02^P))S 

LL = 32sem^(c^(^o2 + ^23/9) - s^(^oi + ^i3/9)) 

Li? = 4m^(-l + /92)((-2 + - 2c0)^i2/9 + (-1 - c,)sec^(^02 + ^23/9) 

+ (l + C0)seS^(^Ol+^13/5))s 

i?0 = 32sem^(c0^o2 - s^Oqi - 0^623/3 + s^9i3f3) 

7. Coefficients of K2 for Mf^^^{LRh'\^+h'\^-) 

Prefactor: g^Qfslys'^/{64tm^) 
00 = se{-C4,eo2 + s^9oi) {f3^ - if (P'' + 2c0 - 3P)s^ 

OL = 2ml,{-l + P^){e^2P{-2{P + If + 2c0{{P + if - s^) + 5^(3 + /3(2 + /?))) 

- S0C^{2sl If + C0{P + lf){0o2 + 02sP) 

+ S0S^{2sl If + C0{(3 + l)')(^oi + 

OR = -2m^(-l + P^)(euP(-2(P - If + 2c0(sl -{(5- if) + s^(3 + (-2 + (5)(5)) 

+ 80C^{{I3 - If - 2sl + C0{I3 - 1)2)(^02 - ^23/5) 

+ S0S^{{(5 - If - 28l + C0{(3 - l)2)(-^oi + OnP))s 
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LO = -2m^(-l + /32)(^i2/3(2(/3 - 1)^ + 2ce((/3 - 1)' - s^) - 5^(3 + (-2 + 

- ses^iiP - If - 2sl + ceiP lf){9oi + ^i3/9) 
+ sec^iiP - If - 2sl + ceiP - lf){eo2 + ^23/9))^ 
LL = 32sgm^,{cg - (3){c^9o2 - 84,601 - €4,623^ + s^^OisP) 
LR = -lQse{-2 + sj - 2c,) (-c^^02 + 84601)771^1 + (3^) 

RO = -2m2^(-l + P'){6,2P{-2{P + 1)^ + 2ce{{f3 + if - s^) + 5^(3 + f3{2 + (3))) 

+ seC4i2sl + If + CeiP + 1)')(^02 - 623P) 
+ S0S4{2sl -{(3 + If + ce{(3 + l)')(-^oi + ^i3/5))s 
RL = 16se(-2 + + 2cg){-C4do2 + s^^oi)rn^(l + 
RR = -32sem^(ce - /3)(-c^(^02 + ^23/9) + s^(^oi + 6M 



8. Coefficients of K2 for Mf^ZsiLRhiY+h]^-) 

Prefactor: (-^;^ - g{;)g^s'^s^/ (128c^tm^(s - m|)) 
00 = se{-C46o2 + s</,^oi) - l)' + 2ce - 2>f3)s' 

OL = 2m^(-l + /32)(^i2/5(-2(/3 + 1)^ + 2c,((/5 + 1)^ - s^) + sl{3 + 13{2 + /3))) 

- seC4{2sl -{13+ If + ce{l3 + 1)^)(^02 + ^23/?) 
+ seS4{2sl - (/3 + 1)' + ce{(3 + l)')(^oi + 

m = -2m^(-l + /3')(ei2/3(-2(/3 - if + 2ce(s^ - (/3 - 1)') + s^(3 + (-2 + (3)(3)) 

+ s,c^((/3 - If - 2sl + ceiP - lf)i6o2 - 623P) 
+ seS4{{(3 - If - 2sl + ce{(3 - lf){-6oi + 6,3(3))s 
LO = -2m^(-l + (3^){6u(3{2{P - 1)' + 2ce((/3 ^ 1)' - - sl{3 + (-2 + (3)(3)) 

- ses^iiP - If - 2s\ + c,(/5 - lf){6ox + ^i3/9) 
+ s,c^((/3 - If - 2sl + ce{(3 - lf){6o2 + 623f3))s 

LL = 32sg7n^{ce - /5)(c^6'o2 - S4601 - C4623P + S4613P) 
LR = -16se{-2 + - 2c,) (-c^^02 + s/oi)^^(l + /9') 

RO = -2m^(-l + /32)(^i2/3(-2(/3 + 1)^ + 2c,((/5 + 1)' - si) + s^(3 + (3{2 + (3))) 

+ seC4{2sl -{P + If + ce{(3 + lf){eo2 - ^23/?) 
+ seS4{2sl + If + ce{(3 + l)')(-^oi + 
RL = 16se(-2 + sl + 2cg){-C46o2 + s^^oi)rn^(l + P^) 
RR = -32sem^{ce - /3) (-0^(^02 + ^23/9) + s^(^oi + ^is/?)) 
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9. Coefficients of K2 for Mf^^^{LRh-^+h-^-) 



Prefactor: g^Qfs'^s'^/{64um^) 

00 = sei-c^9o2 + s^Ooi) {P^ - if {2ce + SP - P')s^ 

OL = 2m^(-l + /3^){euP{-2{P - if + 2ce((/5 - l)^ - s^) + ^^(S + (-2 + f3)P)) 

- SgC^i2sl -iP- I f + CeiP - iy){eo2 + 023^) 

+ ses42sl -iP- 1)2 + ceiP - l)2)(^oi + OisP))s 
OR = -2m^(-l + P^){euP{-2{P + If + 2ce{sl - {P + if) + s^(3 + P{2 + P))) 

+ sec^iiP + If - 2sl + ce{P + lf)%2 - ^23/^) 
+ ses^{{P + If - 2sl + CeiP + l)')(-^oi + enP))s 
LO = -2m^(-l + P^)(euP{2{P + If + 2ce({P + if - s^) - s^(3 + /3(2 + /?))) 

- ses^iiP + If - 2sl + CeiP + l)')(^oi + 0,sP) 
+ sec^iiP + If - 2sl + ceiP + 1)')(^02 + ^23/3))s 

LL = 32sgmfyic0 + /?)(c^6'o2 - s^Oqi - c^923p + s^9i3p) 
LR = -16se(-2 + sl- 2ce)i-c^0o2 + s/oi)mti.(l + P^) 

RO = 2m^(-l + /?2)(^i2/3(2(/3 - 1)^ + 2ceisl - (/? - 1)') - s^(3 + (-2 + P)P)) 

+ sec^i2sl -iP- If + CeiP ' l)')(-^02 + ^23/?) 
+ ses^i2sl -iP- If + c,(/3 - l)2)(^oi - ^i3/?))s 

i?L = lQsei-2 + sl + 2ce)(-c/o2 + sAi)^^(l + Z^") 
= -32sem^(ce + /3) (-0^(^02 + ^23/3) + 5^(^01 + ^is/?)) 



10. Coefficients of K2 for Mf^Zsi^Rh-^+h^^-) 

Prefactor: ig{ + g{^)g^s'^^s^ / (128cLMm,^(s - m|)) 
00= se(-c/o2 + sAi) (/S^ - 2c, - 3/5)^2 

OL = -2m^(-l + /?2)(^i2/3(-2(/3 - 1)' + 2ceiiP - if - s^) + 5^(3 + (-2 + P)P)) 

- sec^i2sl -iP- If + CeiP ~ ^f)%2 + O23P) 
+ ses^i2sl -iP- If + ce(/? - l)')(^^oi + e^3P))s 

OR = 2mlri-l + P^)ieuPi-2iP + If + 2c,(s^ - (/3 + l)^) + ^^(3 + /3(2 + P))) 

+ sec^iiP + 1)2 - 2sl + c,(/3 + lf)ie^2 - 023P) 
+ ses^iiP + If - 2sl + CeiP + l)')(-^oi + OisP))s 
LO = 2m^(-l + P^)i9,2Pi2iP + if + 2c,((/3 + if - s^) - s^S + /3(2 + /?))) 

- ses^iiP + If - 2sl + CeiP + l)')(^oi + O^^P) 
+ sec^((/3 + 1)2 - 2sl + c,(/3 + lf)ieo2 + ^23/9))s 

LL = 32sem^(ce + /3)(-C06'o2 + s<^6'oi + £^6*23/9 - s<^6'i3/5) 
Li? = msei-2 + sl- 2ce)i-c^eo2 + s^^oi)mt^(l + P^) 

RO = 2m^(-l + P^)iei2Pi-2iP - If + 2ceiiP ' ^f ' 4) + ^^(3 + (-2 + P)P)) 

+ sec^i2sl -iP- If + CeiP - ^fM2 - O23P) 
+ ses^i2sl -iP- If + ceiP - l)')(-^^oi + ^i3/9))5 
RL = -16se(-2 + sl + 2ce)i-c^eo2 + s^^oi)m^(l + 
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RR = 32s0m^{ce + /3) (-0^(^02 + O23P) + s^{9oi + 
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